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基于 深度 学 习 的 太阳 活动 区 检测 与 跟踪 方法 研究 
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(1. 昆明 理工 大 学 信息 工程 与 自动 化 学 院 / 云 南 省 计算 机 技术 应 用 重点 实验 室 ， 云 南 昆明 650500; 2. 中 国 科学 院 太 阳 活 
动 重点 实验 室 ， 北 京 100101) 


摘要 : 太阳 活动 区 是 各 类 太阳 活动 的 主要 能 量 来 源 , 剧烈 的 太阳 活动 会 直接 影响 人 类 的 
生存 环境 ， 因此 准确 地 检测 与 跟踪 太阳 活动 区 对 监控 和 预报 空间 天 气 非常 重要 。 本文 基于 深 
度 学 习 框 架 的 YOLOv3-spp 和 DeepSort 提出 了 一 种 太阳 活动 区 检测 和 跟踪 方法 (Active 
Regions Detection and Tracking Method, ARDTM), 该 方法 较 好 地 解决 了 传统 图 像 处 理 方法 吻 
将 一 个 太阳 活动 区 误 检 测 为 多 个 , 或 者 多 个 太阳 活动 区 误 检测 为 一 个 活动 区 的 问题 ; 及 时 捕 
获 到 新 产生 的 太阳 活动 区 和 终止 跟踪 消失 的 太阳 活动 区 ,有 效 地 提高 了 太阳 活动 区 的 跟踪 准 
确 率 。 实 验 结果 表明 该 方法 可 以 较 好 地 检测 和 跟踪 不 同 望远镜 、 不 同时 间 间 隔 序 列 图 像 中 的 
太阳 活动 区 。 
关键 词 : 太阳 活动 区 ; ”检测 与 跟踪 ; ”深度 学 习 ; YOLOv3-spp;  DeepSort 
中 图 分 类 号 : TP391.413 


太阳 活动 区 是 太阳 强 磁场 区 域 ,是 太阳 活动 的 主要 能 量 来 源 。 剧烈 的 太阳 活动 会 导致 恶 
劣 的 空间 和 天气， 比如 瘤 班 和 日 晃 物 质 抛射 等 ， 会 对 地 球 上 电磁 通信 、 电 力 系统 、 无 线 电 传输 
等 产生 不 利 影响 。 因 此 对 太阳 活动 区 进行 检测 与 跟踪 ， 分 析 其 演化 规律 ， 对 人 类 的 空间 探索 
和 地 球 生活 有 着 重要 的 意义 “”。 

目前 用 来 解决 太阳 活动 区 检测 与 跟踪 的 方法 主要 采用 传统 的 图 像 处 理 技 术 ”" 。 比 如 检 
测 主 要 采用 强度 阔 值 法 或 区 域 生 长 等 方法 ， 这些 方法 一 般 需要 设置 参数 ， 比 如 强度 闵 值 、 开 
闭 算 子 阐 值 和 区 域 生长 的 边界 阔 值 。 这 些 参数 对 于 这 些 传统 的 方法 来 说 至 关 重 要 ， 文 [2, 6] 
对 参数 进行 了 详细 讨论 , 并 给 出 了 设置 参数 的 一 些 准则 。 跟踪 太阳 活动 区 则 主要 根据 红 向 较 
差 自 转 定律 来 预测 其 位 置 ， 然 后 根据 欧式 距离 或 特征 的 相关 性 进行 目标 关联 *"。 

这 些 方法 虽然 较 好 地 实现 了 太阳 活动 区 的 检测 与 跟踪 , 但 是 仍 存 在 一 些 问 题 。 比 如 , 文 
[2-4] 的 方法 容易 导致 一 个 太阳 活动 区 被 误 检测 为 多 个 ， 或 者 多 个 相 邻 的 太阳 活动 区 被 误 检 
测 为 一 个 , 从 而 导致 误 跟踪 ; SolarMoni tor ^" 使 用 的 太阳 活动 区 美国 国家 海洋 和 大 气管 理 局 
(National Oceanic and Atmospheric Administration，NOAA) 编号 被 广泛 使 用 ， 但 存在 未 
能 及 时 标注 新 浮现 的 太阳 活动 区 和 仍 标注 了 已 经 在 全 日 面 像 上 消失 的 太阳 活动 区 等 现象 。 

近年 来 ， 随 着 深度 学 习 呈 的 普及 和 CNN 等 神经 网 络 的 发 展 ， 一 些 基 于 深度 学 习 的 目标 检 
测 和 跟踪 的 算法 被 提 了 出 来 .比如 , 典型 目标 检测 算法 有 :R-CNN ^ Faster R-CNN ^, YOLO*, 
YOLOv3 等 。 其 中 ,YOLOv3 借 鉴 残 差 网 络 的 思想 采用 了 新 的 特征 提取 网 络 Darknet-53, 并 结 
合 多 尺度 特征 进行 目标 检测 ， 取 得 了 较 好 的 检测 性 能 。 典 型 多 目标 跟踪 方法 有 : EAMIT 
STAM”, DeepSort 等 。 其 中 ，DeepSort 具 有 较 好 的 跟踪 性 能 和 较 高 的 实时 性 能 。 在 天 文 
领域 的 研究 中 ， 越 来 越 多 的 深度 学 习 方 法 被 采用 。 

本 文 针对 太阳 活动 区 的 演化 特性 , 采用 Y0LOv3-spp 作 为 检测 算法 和 DeepSort 作 为 跟踪 算 
法 并 进行 改进 ， 提 出 了 一 种 太阳 活动 区 检测 和 跟踪 方法 (Active Regions Detection and 


* 基金 项 目 : 国家 自然 科学 基金 ( 11763004，11573012，11803085，U1931107) ; 中 国 科学 院 太 阳 活 动 重点 实验 室 开放 课题 
( KLSA201815) 资助 . 
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Tracking Method，ARDTM) 。 该 算法 可 以 有 效 地 检测 并 跟踪 不 同时 间 间 隔 序列 图 像 中 的 太阳 
活动 区 。 第 一 节 介 绍 了 所 使 用 的 数据 和 制作 的 数据 集 ， 第 二 节 简 述 了 ARDTM 方 法 ， 第 三 节 是 
实验 结果 与 分 析 ， 第 四 节 对 本 文 进行 总 结 与 展望 。 


1 数据 和 数据 集 


本 文采 用 了 太阳 动力 学 天 文 台 (Solar Dynamics Observatory, SDO) 卫星 的 日 震 及 磁场 
成 像 仪 (Helioseismic and Magnetic Imager，HMI) 与 日 球 天 文 台 (Solar and Heliospheric 
Observatory, SOHO) 的 迈克 逊 多 善 勒 成 像 仪 Michelson Doppler Imager, MI) 两 个 设备 的 
太阳 全 日 面 纵向 磁 图 作为 数据 来 源 , 制作 了 太阳 活动 区 检测 与 跟踪 的 深度 学 习 的 训练 和 验证 
数据 集 。 

表 1 列 出 了 分 别 用 于 太阳 活动 区 检测 与 跟踪 的 数据 集 。D1 共 选 取 了 3017 张 全 日 面 图 像 
用 作 制 作 太 阳 活 动 区 检测 网 络 的 训练 和 验证 集 ， 共 选取 了 16 361 个 太阳 活动 区 样本 ， 其 中 
80% 作 为 训练 样本 用 于 训练 深度 学 习 网 络 模型 , 20% 作 为 验证 样本 在 训练 过 程 中 用 于 防止 过 拟 
合 的 验证 。D2 共 选 取 了 33 515 个 太阳 活动 区 样本 用 于 跟踪 过 程 中 重 识别 网 络 的 训练 和 验证 ， 
将 太阳 活动 区 局 部 区 域 缩放 为 128X128 像素 尺寸 ， 同 样 80% 用 于 训练 ，20% 用 于 验证 。 本 文 
根据 SolarMonitor 网 站 上 标注 的 太阳 活动 区 NOAA 编写， 结合 人 工 辅 助 判断 的 方式 ， 采 用 
labelImg 标注 软件 标注 了 数据 集 。 
表 1 太阳 活动 区 检测 和 跟踪 数据 集 

Table 1 The dataset of active regions detection and tracking 

数据 用 途 数据 集 编号 设备 来 源 的 时 间 / 间 隔 样本 数量 


D1 HMI 检测 2011. 01. 01-2017. 12. 31 16 361 


IT RI 


训练 / 验证 
2011. 11. 01-2011. 11. 30/12 分 钟 
集 D2 HMI 跟踪 33 515 
2017. 03. 20-2017. 04. 20/12 分 钟 


2 方法 


ARDTM 算法 主要 分 为 两 个 部 分 : 检测 与 跟踪 。 在 检测 部 分 ， 本 文采 用 了 YOLOv3-spp 模 
型 。 该 模型 的 网 络 结构 如 图 1 所 示 , 由 Darknet-53( 图 1. a) 与 yolo 层 (图 1.b) 两 部 分 组 成 ， 
分 别 用 于 提取 图 像 特征 与 多 尺度 预测 。 方 法 步骤 主要 如 下 : OD) 将 全 日 面 图 像 压 缩 为 608 X 
608 的 尺寸 输入 网 络 。 (2) 输入 图 像 经 过 一 个 卷 积 单 元 (Convolutioan1) 得 到 32 个 608X 608 
的 特征 图 。 卷 积 单元 包含 卷 积 、 批 量 归 一 化 和 leaky relu 激活 函数 三 个 部 分 。(3) 特 征 图 经 
过 步 长 为 2、64 个 卷 积 核 的 卷 积 单元 实现 降 采 样 ， 得 到 64 个 304X 304 的 特征 图 ， 然 后 送 入 
残 差 单 元 (Residual)， 输 出 64 个 304X304 的 特征 图 (1x 残 差 块 ) 1x 表示 此 残 差 块 由 1 个 
残 差 单元 组 成 。 残 差 单元 是 将 输入 与 两 个 卷 积 单元 进行 残 差 操作 。(\4) 上 一 步 得 到 的 特征 轿 
经 过 降 采 样 和 2x 残 差 块 得 到 128 个 152X152 的 特征 图 。 其 中 2x 表示 此 残 差 块 由 2 个 残 差 
单元 组 成 。 (5) 接 下 来 依次 执行 降 采 样 、8x 残 差 块 、 降 采样 、8x 残 差 块 、 降 采样 ，4x 残 差 
块 操作 。 三 个 残 差 块 依次 得 到 256 个 76X76、512 个 38X38 和 1024 个 19X19 的 特征 图 。 
(6) 将 8x、8x 和 Ax 三 个 残 差 块 的 结果 送 入 yolo 层 ( 图 1.b) 。(7)19X19 的 特征 图 经 过 金字 
塔 卷 积 集 (convolutional-spp Set) 等 操作 得 到 第 一 次 的 检测 结果 (0utput1); 同时， 将 金字 
塔 卷 积 集 的 结果 上 采样 至 38X 38 与 第 四 个 残 差 块 的 38X38 的 特征 图 融合 ， 经 过 卷 积 
(Convolutional Set) 等 操作 得 到 第 二 次 的 检测 结果 (Output2); 依次 类 推 ， 得 到 第 三 次 的 检 
测 结果 (Output3) 。 其 中 ， 卷 积 集 操作 包含 5 个 卷 积 单元 ; 金字 塔 卷 积 集 是 在 卷 积 集 基 础 上 
添加 了 空间 金字 塔 池 化 “。 (8) 将 以 上 三 次 的 检测 结果 使 用 非 极 大 抑制 算法 去 除 交 并 比较 大 
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的 预测 框 ， 保 留置 信 得 分 较 高 的 预测 框 作为 太阳 活动 区 的 检测 结果 。 该 网 络 的 损失 值 
部 分 合成 : 目标 定位 偏 移 量 损失 、 目 标 置信 和 度 损 失 和 目标 类 别 损失 。 其 中 ,目标 置信 和 度 损失 
和 类 别 损失 采用 了 二 元 交叉 焙 损 失 ， 目 标定 位 损失 则 采用 了 均 方 误差 进行 计算 。 
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: Type Filters Size Output 

1 i : Convlutional Convolutiona 

i Convolutioana 32 3x3 608x608x32 i Set | 3x3 

i Convo lut ioana 64 3x3/2  304x304x64 n | 4 

i Convolutioana 32 xl 304x304x32 | - Concat Conv : 

i lX|Convolutioanal 64 3x3 304x304x64 | ; x 个 ixl 站 output3| l 

1 Residua | : Up Sampling 

1 Convolutioana 128  3x3/2 152x152x128 : 个 

1 Convolutioana 64 xl 152x152x64 | Convolutional 

| ?x|Convolutioanal 128 3x3 152x152x128 1x1 

1 Residua 1 | ， f 

i Convolutioanal 256 3x3/2 76x76x256 i Convlutional | ,| Convolutiona | 

i Convolutioanal 128 1x1 76x76x128 ME. set 3x3 

! 8x |Convolutioana 206 3x3 76x76x256 | 一 一 i mnm E: 

! Residua | 75 pounds | xl | output2| 

! Convolutioanal 512 3x3/2 38x38x512 Np HPENES QUMMESE — i 

| Convolutioanal 256 1xl 38x38x256 1 | ;LU sumol 1ng i 

! 8x Convolutioana 512 3x3 38x38x512 -一 一 | | 

1 Residua | | Convolutional Eon > Outputl| i 

1 Convolutioana] 1024 3x3/2 19xi9xi024 | B ES i 

; Convolutioana 512 xl 19x19x512 i EACEA nnal 

i 4x s ioana 1024 3x3 19x19x1024 spp Set 3x3 

I I X y 
(a) (b) 


图 1 Yolov3-spp 网 络 结构 示意 图 
Fig. 1 YOLOv3-spp network structure diagram 
TESRRERWUA],. WNAF DeepSort 算法 对 其 进行 改进 ， 主 要 步骤 如 下 : 

CO 根据 太阳 活动 区 的 特性 改进 了 DeepSort 的 特征 提取 网 络 , 本 文 将 其 称 之 为 太阳 活动 
区 重 识别 网 络 ， 如 图 2 所 示 。 首 先 将 检测 到 的 太阳 活动 区 缩放 为 128X128， 经 过 2 次 3X3 
的 卷 积 、1 次 最 大 池 化 以 及 6 个 残 差 块 ， 特 征 图 的 大 小 减 小 到 16X 16, 再 经 过 一 个 全 连接 层 
得 到 128 维 的 全 局 特征 向 量 , 最 后 采用 L2 归 一 化 得 到 特征 向 量 。 在 训练 过 程 中 , 采用 Cosine 
Softmax 函数 得 到 分 类 结果 ， 采 用 交叉 炉 计 算 损失 值 。 需 要 说 明 的 是 检测 部 分 的 网 络 和 大 PB 
活动 区 重 识别 网 络 虽然 都 提取 太阳 活动 区 的 特征 , 但 是 作用 不 同 。 检测 部 分 的 网 络 主要 用 于 
检测 太阳 活动 区 , 而 太阳 活动 区 重 识别 网 络 主要 为 了 判定 太阳 活动 区 在 演化 过 程 中 的 特征 相 
似 性 。 

(2) 用 太阳 活动 区 重 识别 网 络 提取 已 跟踪 到 的 太阳 活动 区 ( 称 为 轨迹 ， 假 设 数量 为 L) 和 
当前 帧 检测 到 的 太阳 活动 区 (假设 数量 为 JD 的 特征 向量 ,计算 之 间 的 余弦 距离 ， 得 到 关联 甜 
WE Cuxi。 然 后 ， 根 据 前 一 帧 所 有 太阳 活动 区 的 位 置 ， 采 用 纬 向 较 差 自转 一” 预 估 在 当前 帧 的 
立 置 。 接 着 ,计算 所 有 预测 位 置 和 所 有 检测 位 置 的 欧 氏 距离 ， 关 过 大 则 将 关联 矩阵 中 对 应 的 
值 设置 为 无 穷 大 。 
(3) ITHIK IE Cxj 采 用 匈牙利 匹配 算法 ”， 得 到 三 类 匹配 结果 :匹配 到 检测 结果 的 轨 
、 未 匹配 到 轨迹 的 检测 结果 和 未 匹配 到 检测 结果 的 轨迹 。 对 于 匹配 到 检测 结果 的 轨迹 则 继 
续 跟 踪 ; 未 匹配 到 轨迹 的 检测 结果 则 为 之 分 配 一 个 新 的 轨迹 ， 当 新 轨迹 连续 3 帧 匹配 到 检测 
结果 时 ， 则 认为 该 太阳 活动 区 产生 ， 否 则 ,删除 该 轨迹 。 对 于 未 匹配 到 检测 结果 的 轨迹 则 提 
取 使 用 纬 向 较 差 自转 预测 结果 的 特征 向 量 ， 与 前 一 帧 对 应 的 特征 向 量 计算 余 粥 距离 。 若 其 值 
在 一 定 范围 内 (本文 实 验 结果 为 0.7) 则 采纳 该 预测 结果 作为 轨迹 的 一 部 分 , 否则 该 活动 区 结 
束 。 但 如 果 连 续 3 帧 均 为 预测 结果 则 也 判定 该 活动 区 结束 。 
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Maxpool 3 
8 (MAX Pooling) 32x64x64 
Patch size:3x3 
Stride:2 
(Convolution) (Convolution) 32x64x64 [M 32x64x64 
Patch size:3x3 Patch size:3x3 Patch size:3x3 
Stride:n-1 Stride:1 Stride:l&l 
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Conv 1 Conv 2 
(Convolution) SR (Convolution) dnd 
Patch size:3x3 Patch size:3x3 

Stride:1 Stride:1 


Residual 4 


Residual 7 
(Residual) 
Patch size:3x3 
Stride:1&l 


Residual 8 
(Residual) 
Patch size:3x3 


Stride:2&l 


64x32x32 64x32x32 


Residual 6 128x16x16 


(Residual) 


Stride:2&l 


Residual 9 128x16x16 Dense 10 iii 


(Residual) L2 
Patch size:3x3 NE normalization 
Stride:1&l 


图 2 太阳 活动 区 重 识别 网 络 网 络 结构 示意 图 


Fig.2 Active regions re-identification network structure diagram 


3 实验 结果 与 分 析 


本 文采 用 HMI 和 MDI 的 四 组 全 日 面 序 列 图 像 进 行 了 测试 。 这 四 组 数据 分 别 代表 了 不 同 的 
观测 设备 和 活动 区 数量 。 表 2 列 出 了 数据 的 基本 信息 , 为 了 便于 描述 , 分 别 用 D3—D6 表示 。 
其 中 ，D3 和 D4 选取 了 MDI 纵向 磁 图 ，D5 和 D6 选取 了 HMI 纵向 磁 图 。D4 和 D6 选取 了 太阳 
活动 区 数量 较 多 的 序列 图 像 , 而 D3 和 D5 选取 了 活动 区 数量 较 少 的 序列 图 像 。 由 于 时 间 间 隔 
较 长 的 数据 太阳 活动 区 的 特征 变化 大 ,跟踪 难度 较 高 , 因此 四 组 数据 均 采 用 时 间 间 隔 较 长 的 
序列 图 像 (为 了 便于 与 SolarMonitor 结果 进行 比较 ， 数 据 与 SolarMonitor 提供 的 数据 时 间 
基本 一 致 ) 。 

K 2 列 出 了 四 组 数据 集 的 跟踪 结果 ， 并 与 NOAA 进行 了 对 比 。 包括 SolarMonitor 上 标注 
的 NOAA 活动 区 数量 、ARDTM 方法 标注 的 太阳 活动 区 数量 、ARDTM 结果 与 NOAA 结果 相同 的 活 
动 区 的 数量 (ARDTMmNOAA) ~ NOAA 活动 区 数量 与 交集 的 差 值 (NOAA- (ARDTM 0 NOAA) ) 、ARDTM 
活动 区 数量 与 交集 的 差 值 (ARDTM- (ARDTMmNOAA) ) 、 交 集 与 NOAA 活动 区 数量 的 百分比 
( (ARDTM N NOAA) /NOAA) 。NOAA 一 共 标 注 了 806 个 活动 区 ，ARDTM 标注 了 872 个 ,两 者 的 交集 
为 755 个 。 ARDTM 方法 标注 的 活动 区 比 NOAA 的 多 一 些 , 主要 由 以 下 四 种 情况 产生 : (1) ARDTM 
及 时 标注 了 新 出 现 的 活动 区 ， 而 NOAA 存在 未 能 及 时 标注 的 现象 ;，(2) ARDTM 较 完 整地 标注 了 
活动 区 的 演化 末期 ， 而 NOAA 已 经 提前 结束 了 对 其 标注 ，(3)ARDTM 标注 了 一 些 没有 被 NOAA 
标注 的 活动 区 ， 这 些 活动 区 的 特点 均 为 磁场 强度 较 弱 、 磁 结构 特征 不 太 明 显 , 但 日 面 已 有 黑 
子 产 生 且 持续 了 3 帧 以 上 ; (4) 一 些 已 具有 磁场 特征 但 未 形成 黑子 的 活动 区 被 ARDTM VR E e 
这 三 种 情况 所 占 比 例 分 别 为 52%、10%、32% 和 6%。 这 意味 着 一 半 以 上 的 情况 是 ARDTM 及 时 标 
注 了 新 出 现 的 活动 区 。 反 过 来 ，NOAA 标注 的 活动 区 比 交 集 多 了 51 个 ， 这 意味 着 ADM 没有 
标注 51 个 被 NOAA 标注 的 活动 区 。 造 成 此 情况 的 主要 原因 有 : (1) 太阳 活动 区 的 特征 已 经 在 纵 
向 磁 图 上 消失 但 仍 被 NOAA 标注 ; (2) ARDTM 漏 标注 一 些 磁 场 强度 较 弱 且 特 征 不 明显 的 活动 区 。 
两 种 情况 所 占 比例 分 别 为 88% 和 12%。 在 这 四 组 数据 中 ，NOAA 在 DA 数据 集 上 及 时 地 跟踪 了 
新 活动 区 和 终止 跟踪 消失 的 活动 区 ， 所 以 本 文 方法 和 NOAA 的 结果 相似 度 最 高 。 


为 了 说 明 以 上 几 种 情况 , K 3 详细 列 出 了 De 测试 集中 前 十 天 的 跟踪 结果 , 第 2—6 列 的 
数据 含义 与 表 2 相同 ， 第 7 列 详细 列 出 了 ARDTM 方法 与 NOAA 标注 结果 不 同 的 具体 情况 。 比 
较 典 型 的 是 AR12185 于 10 月 5 日 出 现 , 但 NOAA 于 7 日 才 开 始 标注 ，AR12183 左 侧 有 个 活动 
区 ,于 10 月 5 日 形成 了 黑子 , 一 直到 14 日 才 消 失 , 但 该 活动 区 一 直 未 被 NOAA 标注 ; AR12173 
已 于 10 月 3 日 消失 于 日 面 ， 但 NOAA 仍 做 了 标注 

表 2 跟踪 结果 
Table 2 The tracking results 


o 


数据 (ARDTM 
. 数据 NOAA ARDTM ARDTM NOAA-(ARDT ARDTM- (ARD 
T IER NOAA) /NO 
来 源 结果 结果 M NOAA MmNOAA) TMN NOAA) 
号 AA 
D3 MDI 2005. 03. 01-2005. 03. 31/24 小 时 68 72 62 6 10 91. 18% 
D4 MDI 2000.07. 12-2000. 08. 11/24 小 时 410 415 398 12 17 97. 07% 
D5 HMI 2010. 11. 05-2010. 12. 04/24 小 时 103 108 91 12 17 88. 35% 
D6 HMI — 2014. 10. 01-2014. 10. 31/24 小 时 225 277 204 21 73 90. 67% 
合计 806 872 755 51 117 91. 82% 


de 3 D6 数据 集 部 分 跟踪 结果 
Table 3 The tracking results of D6 dataset 


h = 3 NOAA ARDTM ARDTM NOAA- (ARDTM ARDTM- (ARDTM 
e 其 备注 


结果 结果 ANOA N NOAA) N NOAA) 
C 
nie ARI2171 已 从 全 日 面 消失 ，NOAA 仍 标记 ;NOAA 
p 2014-10-1 10 10 9 1 
QJ) 未 及 时 标注 AR12182; 
LER AR12175 已 从 全 日 面 消失 ，NOAA 仍 标记 ; NOAA 
a» 2014-10-2 9 9 8 1 
EN 未 及 时 标注 AR12182; 
AR12173 已 从 全 日 面 消失 ，NOAA 仍 标记 ; NOAA 
2014-10-3 8 10 7 3 
未 及 时 标注 ARI2182. ARI2183. ARI2184; 
g 2014-10-4 10 9 9 0 ARI2172 已 从 全 日 面 消失 ，NOAA 仍 标记 ; 
E: = AR12176 已 从 全 日 面 消失 ，NOAA 仍 标 记 ; NOAA 
2014-10-5 9 10 8 2 未 及 时 标注 AR12185; ARDTM 标注 但 NOAA 未 标 
iE 1 个 活动 区 
AR12176 和 AR12180 已 从 全 日 面 消失 ，NOAA 仍 
2014-10-6 9 9 7 2 2 标记 ; NOAA 未 及 时 标注 AR12185; ARDTM 标注 但 
NOAA 未 标注 1 个 活动 区 ; 
NOAA 未 及 时 标注 AR12186; ARDTM 标注 但 NOAA 
2014-10-7 8 10 8 0 2 E 加 
未 标注 1 个 活动 区 ; 
2014-10-8 9 10 9 0 1 ARDTM 标注 但 NOAA 未 标注 1 个 活动 区 
AR12181 已 从 全 日 面 消失 ，NOAA 仍 标记 ; ARDTM 
2014-10-9 8 8 7 1 1 
标注 但 NOAA 未 标注 1 个 活动 区 
AR12178 和 AR12179 已 从 全 日 面 消失 ，NOAA 仍 
2014-10-10 6 5 2 1 g mu " £ 
标记 ; ARDTM 标注 但 NOAA 未 标注 1 个 活动 区 
合计 87 91 71 10 14 


ARDTM 较 好 地 解决 了 文 [4 中 提 到 的 一 个 太阳 活动 区 被 误 检测 为 多 个 ， 或 者 多 个 太阳 活 
动 区 被 误 检测 为 一 个 活动 区 的 问题 。 图 3 (8) RH T CLA P 11(a) 的 MWDI 同时 间 数 据 (2002 
FE5 月 4 日 12:480D ,红色 实 线 框 为 ARDTM 的 检测 结果 , 蓝 色 虚线 框 为 文 [4 的 结果 。 其 中 ， 
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1 和 2 两 个 区 域 分 别 为 AR09934 和 AR09936， 而 文 [ 币 误 把 他 们 识别 成 了 一 个 活动 区 。(b) 为 


HMI 的 2010 年 11 H 12 日 数据 。 同样 ，3 和 4 分 别 为 AR11123 和 AR11121， 而 用 文 [4] 方 法 


误 识别 成 了 一 个 活动 区 。(c) 为 HMI 的 2010 年 11 月 7 日 的 AR11120， 
活动 区 ， 但 采用 文 [种 方法 将 其 检测 成 了 两 个 活动 区 。 


x 


ARDTM 顺利 检测 为 一 个 


2002/05/04 12:48:00UT 2010/11/12 19:12:00UT 2010/11/07 23:12:00UT 


图 3 ARDTM 的 检测 结果 图 。 


Fig.3 The detection results images diagram of ARDTM 


图 4 展示 了 D6 数 据 集中 连续 三 帧 的 跟踪 结果 (a) 、(b) RI (c), 时 间 分 别 是 2014 年 10 月 20 日 


22:48UT、21 日 23:00UT 和 22 日 22:48UT, 图 中 粗 框 代表 了 检测 到 的 太阳 


活动 区 匹配 到 了 轨迹 ; 


虚 框 表示 该 区 域 未 被 检测 或 未 被 跟踪 匹配 , 而 采用 预测 机 制 填补 的 结果 。 跟踪 框 的 左上 和 角 数 


字 表 示 ARDTM 对 太阳 活动 区 的 跟踪 编写。 通过 与 对 应 时 间 的 NOAA 对 比 ， 


ARDTM 具 有 良好 的 跟踪 


结果 。 比 如 ，(1) ARDTM 及 时 捕获 到 了 新 的 活动 区 。ARDTM 在 (a) 及 时 


B3k T ARI2194, TjNOAA 


在 (c) 才 被 标注 ，ARDTM 比 NOAA 提 前 了 2 天 。NOAA 中 AR12183、AR12186 


、AR12189、AR12190、 


AR12191、AR12192 和 AR12195 这 几 个 活动 区 ，ARDTM 均 比 NOAA 更 早 检 测 到 并 及 时 开始 跟踪 。 
(2)ARDTM 及 时 终止 跟踪 消失 的 活动 区 。ARDTM 在 22 日 (c) 及 时 终止 跟踪 AR12189， 而 NOAA 直 到 


26 日 才 终止 跟踪 该 活动 区 。 同 样 ，(a) 图 的 AR12186 也 是 类 似 的 情况 。 


(3) 预测 机 制 较 好 地 弥 


补 了 检测 算法 漏 检 带 来 的 漏 跟踪 或 误 跟踪 问题 。ARDTM 在 Ca) AI (b) 都 通过 预测 机 制 连 续 跟踪 


了 AR12189。(4)ARDTM 检 测 并 跟踪 到 了 NOAA 漏 检 的 活动 区 。(a) 图 中 编 


号 为 23 和 29 的 活动 区 和 


(b) 中 编号 为 33 的 活动 区 ， 在 演化 过 程 中 都 形成 了 黑子 ， 生 命 期 分 别 为 8、6 和 6 帧 ,但 都 未 被 


NOAA 标 注 。 


"ARI2190 


(a) 
图 4 ”ARDTM 在 D6 数据 集中 连续 3 帧 的 跟踪 结果 。(a) 2014 年 10 月 20 H 22: 48UT; (b) 2014 f& 


E10 月 21 H 23: 00UT; (c) 2014 


年 10 H 22 H 22: 48UT。 


Fig.4 ARDTM tracking results of three consecutive frames in the D6 dataset. (a) Oct. 


2014 23:00UT; (c) Oct.22, 2014 22:48UT. 


20, 2014 22:48UT; (b) Oct. 21, 


5 展示 了 图 4 中 AR12194 以 及 D4 数据 集中 ARO9093 两 个 活动 区 的 完整 跟踪 结果 。(a) 


到 (1) 子 图 是 AR12194 从 出 现 到 消失 的 跟踪 过 程 ， 时 间 从 2014 年 10 


H 20 日 22:48:00UT 到 


31 日 19:48:00UT。 该 活动 区 穿越 了 整个 日 面 ， 磁 场 结构 和 面积 变化 都 非常 大 。(n) 到 (t) 子 


图 是 AR09093 从 出 现 到 消失 的 跟踪 结果 ， 时 间 从 2000 年 7 月 16 日 1 


2:48:00UT 到 22 H 


201912.00028v1 


chinaXiv 


ChinaXiv& fERBTIJ 


12:48:00UT. (m) I (u) 是 ARO9093 产生 的 前 一 帧 和 消失 的 后 一 帧 对 应 位 置 图 。 该 活动 区 演化 
周期 时 间 较 短 ， 特 征 变化 也 很 大 。ARDTM 及 时 捕获 且 及 时 停止 了 跟踪 该 活动 区 。 其 中 ，(t) 
时 刻 通 过 预测 机 制 填补 了 漏 检 ，(u) 时刻 及 时 停止 了 跟踪 。 图 6 显示 了 图 5 中 AR09093 和 

AR12194 两 个 活动 区 的 磁 通 量 和 面积 的 变化 曲线 ( 均 已 进行 了 投影 改正 ), 可 以 看 出 其 磁 通 量 
和 面积 都 具有 较 强 的 渐变 性 。 


AR12194 3 AR12194 4 AR12194 5 AR12194 6 , 


AR12194 1 ART12194 2 


x K * + 者 yë n 


014/10/20 2 2480001 014/1021 200 000 2014/10/22 22:48:00UT 2014/10/23 22:48:00UT 2014/10/24 .22:48:00UT 
(a) (b) (c) (d) (e) (f) 
AR12194 7 S AR12194 8 AR12194 9 AR12194 10 AR12194 11 
所 "s B 
' fae E. PP EV P 
æ Ei ^ 5 i 


w A 


ARI2IS4 12 


Í 
2014/10/26 19:48:00UT 个 / 2014/10/28 20;00:00UT 2014/10/29 19:48:00UT 2014/10/30 19:48:00UT 2014/10/31 19 48:00UT 


RS, AR09093 3 AR09093 4 PAS RA 
AR09093 1 z 
* og; 
Ly a 
r jà 
Lo -= » 


2000/07/15 12:48:00UT 2000/07/16 12:48:00UT 2000/07/17 12:48:00UT 2000/07/18 12:48:00UT [2000/07/19 12:48:00UT 2000/07/20 12:48:00UT 


AR09093 6 ARO09093 7 


d 


2000/07/21 12:48:00UT 2000/07/22 12:48:00UT 2000/07/23 12:48 00UT 
(s) (t) (u) 


图 5 (a) 到 (k) 是 D6 数据 集中 AR12194 完整 的 跟踪 结果 。(n) SU (t) 是 D4 数据 集中 AR09093 完整 的 跟踪 结果 。 


Fig.5 From (a) to (k) is the complete tracking result of the number AR12194 in D6 dataset. From (n) to (t) is the 
complete tracking result of the number AR09093 in D4 dataset. 
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(a) 
图 6 (a) 图 5 中 AR12194 磁 通 量变 化 ，(b) 图 5 中 AR09093 磁 通 量变 化 


Fig. 6 (a) The change of magnetic flux of the AR12194; (b) The change of magnetic flux of the AR09093. 


4 总 结 与 展望 


本 文 基于 深度 学 习 提出 了 一 种 太阳 活动 区 检测 与 跟踪 的 方法 ARDTM。 该 方法 首先 使 用 
YOLOv3-spp 网 络 检测 太阳 活动 区 ， 然 后 结合 较 差 自转 改进 DeepSort 算法 跟 踊 太阳 活动 区 。 
经 过 网 络 训练 后 ， 本 文 分 别 用 四 组 数据 进行 了 测试 ， 并 与 NOAA 的 结果 进行 了 比较 ，ARDTM 
跟踪 标注 的 结果 与 NOAA 的 交集 约 为 NOAA 的 92%。 

本 文 的 训练 集 只 采用 了 HMI 数据 ,但 无 论 是 测试 HMI 的 数据 ,还 是 MDI 的 数据 ， 都 能 很 
好 地 检测 和 跟踪 太阳 活动 区 ， 表 现 出 了 较 好 的 泛 化 性 。 为 了 与 SolarMonitor 的 结果 进行 逐 
一 比较 ， 本 文 主要 采用 与 SolarMonitor 网 站 公布 的 数据 时 间 基 本 吻合 的 数据 进行 了 测试 ， 


时 间 间 隔 约 为 24 小 时 。 本 文 还 对 不 同时 间 间 隔 的 HMI 序列 图 像 进行 了 测试 ， 时 间 间 隔 选 取 
为 12 分 钟 、1 小 时 、4 小 时 ……, 一 直到 48 小时。 总 体 来 说 ,间隔 为 12 分 钟 的 数据 由 于 活 
动 区 特征 变化 不 大 ， 跟 踪 结 果 最 好 ; 一 直到 24 小 时 的 时 间 间 隔 跟 踊 效 果 都 还 不 错 ;， 超过 24 
小 时 的 时 间 间 隔 就 会 出 现 较 多 的 误 跟 踪 现 象 。 总 的 来 说 , 该 方法 较 好 地 解决 了 传统 方法 中 一 
个 太阳 活动 区 被 误 检测 为 多 个 ， 或 者 多 个 相 邻 的 太阳 活动 区 被 误 检测 为 一 个 的 问题 。 而 且 ， 
该 方法 能 及 时 捕获 到 新 出 现 的 活动 区 , 亦 能 及 时 终止 跟踪 消失 的 活动 区 , 并 且 通 过 预测 机 制 
弥补 检测 算法 漏 检 带 来 的 误 跟踪 或 漏 跟 踪 问 题 。 当 然 ， 该 方法 也 还 存在 一 定 的 问题 ， 比 如 | 
于 单一 地 根据 磁场 结构 的 特征 检测 和 跟踪 活动 区 , 会 导致 一 些 具有 做 场 结构 但 未 形成 黑子 的 
特征 被 误 跟踪 情况 。 在 未 来 的 工作 中 , 拟 考虑 结合 其 他 波段 观测 的 图 像 进一步 提高 算法 的 疹 


Hr 


确 率 。 
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Abstract: Active regions(ARs) are the primary source of energy for various solar 


activities. The violent solar activities adversely affect human living environment. 


Therefore, accurate detection and tracking of ARs are very important for monitoring 


and forecasting the space weather. In this paper, we propose an AR detection and 


tracking method (ARDTM) based on the deep learning model comprising of the 
YOLOv3-spp and DeepSort. This paper solves the problem that one AR is 
mis-detected as multiple ARs, or multiple ARs are mis-detected as an AR. Besides 
that, it captures the new ARs and terminates the disappeared ARs in time. The method 
improves the precision and recall of detecting and tracking ARs. It can be used for 
detecting and tracking ARs in the solar full-disk longitudinal magnetograms observed 
from different telescopes, or different time interval series-images. 
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